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Testing the accuracy of the overlap criterion 
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Here we investigate the accuracy of the overlap criterion when applied to a simple near-integrable 
model in both its 2D and 3D version. To this end, we consider respectively, two and three quartic 
oscilators as the unperturbed system, and couple the degrees of freedom by a cubic, non-integrable 
perturbation. For both systems we compute the unperturbed resonances up to order 0(e 2 ), and 
model each resonance by means of the pendulum approximation in order to estimate the theoretical 
critical value of the perturbation parameter for a global transition to chaos. We perform several 
surface of sections for the bidimensional case to derive an empirical value to be compared to our 
theoretical estimation, being both in good agreement. Also for the 3D case a numerical estimate is 
attained that we observe matches the critical value resulting from theoretical means. This confirms 
once again that reckoning resonances up to 0(e 2 ) suffices in order the overlap criterion to work out. 

Keywords: Chaos - Resonances - Theoretical and Numerical Methods 



I. INTRODUCTION 



Though the stability problem of Hamiltonian systems has been almost completely elucidated by a rigourous sequence 
of theorems that build up the so called KAM theory (see for instance [l| and [2] , together with the original references 
therein: Q, [J] and |5|), the application of the results of the KAM theory to a specific system is far from being an 
easy task. In fact, it turns out to be much simpler to take advantage of the heuristic Overlap Criterion, which seems 
to provide similar estimations to those resulting from the KAM theory. 

The overlap criterion due to Chirikov (see [l|) has been largely used in many different fields, its probably most 
popular application being the study of inestabilities in the Solar System as well as in other planetary models (see 
for example [7j, |8|, [9(). In any case, since the widespread model for a resonance is the pendulum approximation, 
the overlap criterion applies directly to the intersection of their associated unperturbed separatrices (or heteroclinic 



l2^ intersections). 



iQ , In his pioneer work on the standard map 1] Chirikov shows that the application of the overlap criterion to primary 

£»o ■ resonances overestimates the actual value of the critical parameter, K c , and only when high order resonances are 
considered, does the overlap criterion succeed in providing a more accurate value for K c . In fact, the author shows 
that on including the third harmonics resonances, the overlap criterion leads to K c sa 1, rather close to the empirical 
value. 

In the present effort we address a similar analysis to that performed by Chirikov, but using a 2D and a 3D near- 
integrable Hamiltonian systems, namely, two and three uncouppled quartic oscillators perturbed by a cubic term. The 
. ^h ', 3D version of this model has been studied in [T(| , [12| , where the authors numerically investigate the global dynamical 
S^ ■ properties of the model and estimate the critical value e c beyond which the system is globally chaotic, i.e. for which 
less than the 10% of the energy surface corresponds to invariant tori. By means of the overlap criterion we derive 
such a critical value for the perturbative parameter e on considering not only primary but also high order resonances, 
and the perturbation Fourier series truncated at 0(23~ 2 ) in their coefficients. We then compare, for each case, the 
theoretical critical value with that obtained by numerical means. 

The paper is organized as follows. The dynamical system under study in its 2D version is described in Section 
II and its relevant resonances at 0(e) are obtained in Section III, their widths being determined in Section IV. The 
resonances at 0(e 2 ) are provided in Section V, where an estimate of the critical value of the perturbative parameter 
is provided. For the sake of comparison, an empirical estimate of such a value is obtained in Section VI by recourse of 
performing several surfaces of section for the system. The 3D model is addressed in Section VII, whose resonances at 
order 0(e) and 0(e 2 ) are given in Sections VIII and IX respectively. Section X is devoted to the theoretical estimate 
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of the critical parameter, which is shown to be in good agreement with the one given in [12j |. A final discussion is 
provided in Section XL 

II. THE 2D DYNAMICAL SYSTEM 

Here we will be concerned first with a 2D perturbed quartic oscillator. In cartesian coordinates the system is 
described by the following Hamiltonian (see [lfjl|): 

H(p, q) = \(PI+ Pi) + \^+ V') + &?* (!) 

where e is a perturbative parameter that controls the strength of the perturbation. On setting e = we recover the 
integrable quartic oscillator's Hamiltonian (see [0|,[l3| and references therein), whose solutions are given by 

x(t) = x {h{) J2™=i a n cos (( 2ri - l)oJi(hi)t), 

(2) 
y(t) = Vo(h 2 ) J2^=i a nCos((2n- l)ui 2 (h 2 )t), 

where we have used the following definitions: 

Xo(hi) = 4/3/ii 1/4 , 
y (h 2 ) = 4/% 2 1 / 4 , 

u>i(hi) = V2/3hi 1 ' 4 , i = l,2 ^ 

1 



" cosh ((» - 1/2)tt) ' 

/3 = tt/2K(1/V2) « 0.847, 

where K(k) denotes the complete elliptic integral, and the coeficients in the Fourier expansions ^ satisfy: 

Q-n+i _ J_ 
a n 23 

The third equation in ([3]) reveals the dependence of the frequencies on the unperturbed energies, enabling us to 
get the functional relationship between the latter and the unperturbed action variables, namely, hi = Alt ' , with 
A=(3/3/2 V / 2) 4/3 . 

With this relation in mind, and taking into account that the angle variables are 9i = u)i(hi)t, i = 1,2, the complete 
Hamiltonian, in terms of the action-angle variables of the unperturbed Hamiltonian, can be recast as: 

H(I,9)=Ho(I)+eV(I,e), (4) 

where 

H (I) = AihW + I? 3 ), 
V{I,8) = V(I) J2 ocn m k{cos(2(n + m-l)0i±(2k-l)0 2 )+cos(2(n-m)di±(2k-l)e 2 )} (5) 

n,m,k=l 

with a nm k = a n a m ak, and V(I) = 2 5 / 2 3/? 4 / 1 I 2 , the ± sign meaning that both terms are included in the series. 



III. RESONANCES AT 0(e) 



A glance at the perturbation series given in equation ([5]) reveals that the number of resonant terms at first order 
in the perturbative parameter is unbounded, which is a drawback to take into account the width of every resonance 
at such an order. 

However, the strong dependence of the Fourier amplitudes on (n + m + k), through the quantities a nm k ~ 
l/23( n+m+k ~ 3 \ gives us a good hint on how to gather the 0(e) resonances and where to cut the series, our at- 
tempt being to keep terms only up to 0(l/23 2 ). 

All the possible combinations of n, m and k verifying that n + m + k < 5, yield 24 different vectors which are listed 
in Table [H together with the number of times they appear in a term with coefficient a nm k of a given order in 1/23. 
Thus, Nq denotes the number of times the vector appears with coefficient a nm k — ct\, N\ the number of times it arises 
with coefficient a\a.2(~ a^/23), and N2 corresponds to either the coefficient a\a^ or a\aj, (which are approximately 
a\ /23 2 ). From now on, these vectors will be appointed as harmonics at order 0(e, 1/23 2 ). 
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3 


(-4,1) 








1 


(-4,-1) 








1 



TABLE I: Harmonics in the Fourier expansion © at 0(e, 1/23 2 ). 

On applying the resonance condition m ■ u> = 0, with m E Z 2 /{0}, to the unperturbed system, the following 
relation between the energies in each degree of freedom is obtained: 



ra\h 1 



1/4 



/1/4 

m2h 2 



0, 



(6) 



which implies that the resonance structure in energy and action space consists of straight lines (with positive slope) 
given by 






(7) 



and 



T 1 ' — 
1 2 — 



(8) 



respectively. 

Moreover, eq. ([6]) indicates that m\m2 < 0, whence, those vectors having both components with the same sign 
must be discarded. Let us notice however, that not all of the remaining harmonics at order 0(e, 1/23 2 ) are actually 
resonant, as it will be discussed in the forthcoming section. 



IV. WIDTH OF THE RESONANCES AT 0(e) 



On computing the widths of resonances, the pendulum's approximation (see for instance [l|, [l3j) provides a suitable 
description whenever each resonance is assumed to be isolated from the rest. 



Notice should be taken that, before proceeding to estimate the width of a given resonance, all the coefficients a nm k 
associated to the same trigonometric function are to be added together into a single one. Indeed, for each given vector 
m, we should define the coefficient: 

m+n+/c<5 

where n, to, and k are natural numbers that combine to form the vector m in any of the four ways displayed in eq. 
©. Let us remark that a m w a?(7V + iVi/23 + N 2 /23 2 ). 

Inasmuch the pendulum's approximation has been applied to this single resonant term, the new (resonant) Hamil- 
tonian turns out to be: 

H r (pi,ih) = j^+eV{r)a m caaih, (9) 

with 

M" 1 = mj-^-m, = ml-^-, I = I r + m Pl , *=m^, (10) 

where the sum over repeated indexes should be understood. 

Let p r be the maximum variation of p\ within the oscillation regime, then 



p r = 2 (e MV{I r )a nmk ) = 2 7 ' 2 fi 1 l 



1/2 



V2 1/2. an5/8 



e ' a 



urn 



(Kf /S - (ii) 



As a consequence of the simple pendulum dynamics, the maximum displacement of the unperturbed action variables 
depends on m and p r in the fashion: (AJ) r = ( I — I r ) max = p r m. 

Furthermore, the maximum displacement of the unperturbed energy is given by |(Aft,) | = |wj"(AJi)^j|; and it is 
the maximum amplitude attained in the oscillation of any of the unperturbed energies that measures the width of the 
resonance. 

Let us now recall that in any 2D problem the energy conservation condition h = hi + fi2, together with the resonance 
condition m\h-y + TO2^2 = allow both h\ and h T 2 to be written in terms of the total unperturbed energy h. Thus, 
the amplitude can be recast in terms of mi, rrn, e, a m and h as follows: 

|TO?+TO*| n /8 . 12 . 

|(Afta)^| = |(A^)^|. 

The last identity in (fT!?)) is due to the fact that in presence of a single resonance, the motion of the system is tangent 
to the unperturbed energy surface, and for this particular model such a surface is given by h — hi + h 2l which leads 
to A/i 2 = -Afti. 

We note that the width of the resonances at C(e, 1/23 2 ) depends on the harmonic numbers in the manner shown 
in Fig. [H so that those resonances with values of 7 = m<z/mi out of the range [0.5, 2.5] should be narrow. 

A glance at both eq. (fT2)l and Fig. Q] reveals that those harmonic vectors with any of its components equal to zero 
do not change the energies, and consequently, they should not be considered as resonant vectors. 

Further, since the perturbation terms are even, if m is a resonant vector then —m is also a resonant one (both 
corresponding to the same resonance). Hence, for each resonance just one representative resonant vector can be 
considered, encompassing into its concomitant coefficient the contribution corresponding to its opposite vector as 
well. All the relevant data required to compute the width of each resonance at 0(e, 1/23 2 ) is displayed in Table HT1 
where we have also included the value of h\ corresponding to a total unperturbed energy of h = l/(4/3 4 ) jw 0.485 (the 
one used in |10j). 

We have computed the resonance widths corresponding to e in the range [0, 0.5]. Fig. [2]displays both the maximum 
and minimum values of hi for each resonance vs. the perturbative parameter e; the total unperturbed energy being 
h = l/(4/3 4 ). We observe that for e ~ 0.15, the (6, — 1), (4, — 1), (2 — 1) and (4,-3) resonances do overlap, but lie 
far away from the (2, —3) and (2, —5) resonances. Therefrom we could infer that the energy surface presents two 




FIG. 1: Dependence of the resonance width on m2Jm\. 
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0.0019 


(6,-1) 
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0.0004 



TABLE II: Resonant vectors at 0(e) up to 0(l/23 2 




FIG. 2: Widths of those resonances at 0(e) up to 0(1/23 ), in terms of e. 



unconnected regions of chaotic motion, so that a global transition to chaos does not take place for e < 0.5, which 
leads to a critical thcorical value for the perturbation parameter e c 3> 0.5. 



V. RESONANCES AT C(e 2 ) 

In those regions of phase space which are far from any primary resonance (i.e. where the diophantic condition 
holds for every primary resonance), we can introduce new cononical variables, (J,ip), in such a fashion that the 
transformed Hamiltonian consists of a part depending on the new momentum and a perturbation that, though being 
non-integrable, has an amplitude of 0(e 2 ). 

2 (id a H (j)a9(j,(p)a<t(j ) <p) , dv(j,cp)d<s>(j,<p) 



nj, V )=H (j) + <>{- dJ £ Jt > — — + — d ;r ) + ° {e% (13) 

where $ stands for the trigonometric part of the generatrix function of the canonical transformation: 
F(J,d) = J-6 + e$(J,d) $(J,0)= Y^ $m(J)sin(m-0), 

m€Z 2 /{0} 

with $ m (J) = -V m (J)/m ■ w(J). 

After computing the right side of expression (|13[) one finds: 

W(J, <p) = H (J) + e 2 J2 C(m ' m '' J ) i cos (( m + m ') ' v ) + cos (( m ~ m "> ' V)} + °^ 3 ^ ( 14 ) 

where the coefficients are given by 

{/ , T 2/3 T 2/3 . , T 4/3x / n , T l/3 7 2/3 . , 7 4/3 t-1/3n ^ 
o3o4/3 Q 28/3 ( m l TO l-V -V +m 2 771 2 J 1 / ) 4o Q8 ( 2m l J l J 2 + m 2 -V ^2 ) I 
1 A p — ? T7TT7 — ; — 7~rr\ 2 6 p ? — ; — nv\ r ■ 
(m ■ uj(J)){m' ■ <jj(J)) (m' ■ u){J)) J 

There are several relevant facts to be remarked: (i) on working at 0(e) we have only considered resonances up 
to 0(l/23 2 ) in the Fourier coefficients; therefore, the series in equation (|14[) should actually be replaced by a finite 
sum over those harmonics, m and m', whose associated coefficients {ot m and a m >, respectively) are such that their 
product is of order either a\ , af/23 or af/23 2 ; (ii) there are many different pairs of harmonics (m, m') at 0(e), which 
combine into the same harmonic n at 0(e 2 ); (iii) as a consequence of the evenness of the perturbation term, if n is a 
resonant vector — n is also a resonant one; (iv) the resonance condition implies that n\fi<i < 0; and (v) the condition 
of being far from resonances at 0(e) implies that we must discard all those harmonics which are a multiple of any 
resonant vector at 0(e, 1/23 2 ). 

To cope with the situation set up by the issues (ii) and (iii), wc have added all the concomitant contributions into 
a single coefficient D, namely, 

D(n,J) = Y, C(m,m',J), (15) 

m,m' 

where the sum extents to all the harmonics (m, m') at 0(e, 1/23 2 ) such that n = m+m', n = mm', —n = m+m', 
or — n = m — m! , and for which 0(a m a m >) is not greater than a^/23 2 . 

Taking into account all the above mentioned considerations, the Hamiltonian can be written in the form: 

H(J, if) = H (J) + e 2 J2 D ( n ' J ) cos ( n ' V)- (16) 

Therefore, in the vicinity of a resonant torus J r , and by recourse of the pendulum approximation, we obtain the 
new resonant Hamiltonian: 

« r (Pi,*i) = ?-+«ocos(tti), (17) 

2(1 

where 



U {) 



duV 

e 2 B(n,J r ) 



Thus, the maximum displacement of the unperturbed action variables is given by (A J) r = ( J — J r ) r 
with V r = 2 (/x|Wo|) , and the widths of the resonances at 0(e 2 ) are given by 



(18) 
V r n, 



iw;hW \| , ; l L 

(n-J + rig) ' 



Z a ^ 



2(2|ni| 3 m' 1 + |ri2| m' 2 ) n\m! 1 m\ + n\m! 2 m2 



± 



|niTl2|(m. / 1 |7l2| + m 2 |ni|) ( TO ll n 2| + m' 2 |ni|) 



1/2 



where the plus sign corresponds to n = m + m' and the minus sign to n = m — m! , m, ought to be written in terms 
of m'i and n,. As a consequence of this last expression, resonant vectors are compelled to have no null components. 



vector 
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(2,-6) 

(2,-4) 
(2,-2) 
(6,-4) 
(6,-2) 


0.4791 
0.4565 
0.2425 
0.0800 
0.0059 



TABLE III: Resonant vectors of C(e 2 ) up to C(l/23 2 




FIG. 3: Width of the resonances up to 0(e 2 , 1/23 2 ) vs. e. 

The harmonics satisfying all the stated conditions for the perturbation at order e 2 turn out to be just five, which are 
listed in Table lim They will be referred to as resonant vectors at 0(e 2 , 1/23 2 ). We have computed the concomitant 
resonance widths for e in the range [0, 0.5] . The results are presented in Fig. G2 where also the resonances corresponding 
to 0(e) have been included. Let us recall that the adopted value for the total unperturbed energy is h = 1/4/3 4 . 

From Fig. [3] and Table ULTl we notice that the arising of the (2, —2) resonance connects the two sets of resonances 
that at C(e, 1/23 2 ) appeared isolated for the considered e range. The remaining resonances at 0(e 2 , 1/23 2 ) appear 
completely overlapped with either set of resonances at 0(e, 1/23 2 ). From this plot, we could derive the critical value 
for the perturbative parameter, e c ~ 0.28. 



VI. NUMERICAL ESTIMATION OF THE CRITICAL VALUE OF THE PERTURBATION 

PARAMETER 

In this section we empirically estimate the value of e c by means of Poincare Surfaces of Section (SOS). To this aim, 
we take the intersections on the plane x = (actually |x| < 10 -8 ) whenever p x > 0, for several initial conditions along 
the y-axis. 

Fig [4] displays the SOS's corresponding to e = 0.12 (on the left) and to e = 0.14 (on the right), respectively. There 
we can distinguish the (2, —1) resonance, very close to the last invariant curve that corresponds to the y-axis periodic 
orbit - (1,0) resonance-, and the (2, —2), (2, —3) resonances as well. The (2,-4), (2,-5) and (2,-6) resonances do 
not show up due to the fact that they are completely distroyed by overlap, as could be seen from Fig. [3] for this value 
of e. It is important to remark that several higher order resonances do appear which have not been theoretically 
computed. 

For e = 0.12, we observe that the chaotic domain where the (2, —4), (2, —5) and (2, —6) resonances appear distroyed 
by overlap is bounded by some KAM tori and thus, it remains unconnected with the outer chaotic component around 
the (2,-3) and (2,-1) resonances. On the other hand, for e = 0.14, both chaotic zones are connected, leading to a 
global transition to chaos, in the sense that any orbit could explore almost all the chaotic component of phase space. 
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FIG. 4: Poincare surfaces of section for e = 0.12 (on the left), and for e = 0.14 (on the right). 

Therefore, it looks like e c lies somewhere in the range (0.12,0.14). After performing a rather thorough numerical 
exploration, we have noted that for e = 0.135, several KAM tori do persit, which are shown in FigEl where a zoom 
in the window [0,0.4] x [0,0.3] is presented. There, such KAM tori can be clearly distinguished and are seen to 
definitively separate both chaotic domains in phase space. Nevertheless, this bounded region of chaotic motion, does 
not involve the resonances we are taking into account to derive the analitical esimation of e c , but high order ones. 

Therefore, from experimental means, we may state that e > 0.135 is a good lower bound for the critical value of 
the perturbation parameter. 

VII. THE 3D MODEL 

Now we will focus on a 3D version of the dynamical system under study. In cartesian coordinates, its Hamiltonian 
is given by: 



H(p, q) = \(p 2 x + pI+P 2 z ) + \{x 4 + y 4 + z 4 ) + ex 2 (y + z). (19) 

Let us notice that for a null value of the perturbative parameter we recover the three independent one dimensional 
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FIG. 5: Poincare surface of section for e = 0.135 illustrating the existence of KAM tori that separate both chaotic components, 
quartic oscillators, whose solutions x(t) and y(t) are the ones given by equation ([2]), while z(t) allows for the expression: 



z(t) = z (h 3 ) ^2 a « cos ((2n - l)iv 3 (h 3 )t) , 



(20) 



n=l 



where z (h 3 ) = 4/3/i 3 1/4 . 

In terms of the action-angle variables of the unperturbed Hamiltonian, the complete Hamiltonian Q19p can be recast 
as: 



with 



H(I,0) = H o (I)+eV(I,9), 

Ho(/)=^(/i 4/3 + / 2 4/3 + / 3 4/3 ), 

OO 

V{1, 6) =V 12 (I) Yl anmk ( cos ( 2 ( n + m - l )°^ ± ( 2k ~ 1 ) 6 ' 2 ) + cos ( 2 ( n _ m ) 6 ' 1 ± ( 2fc ~ 1 ) 6 ' 2 ) } + 

n : m,k— 1 

CSO 

+ ft 3 (I) ^ a nmk {cos (2(n + m - 1)01 ± (2k - 1)6 3 ) + cos (2(n - m)0i ± (2fc - 1)0 3 ) } , 

n.m.k— 1 



(21) 



(22) 



where the new quantities Vij(J) = 2 5 / 2 3/3 4 / 1 L have been introduced; the ± sign meaning that both terms are 
included in the series. 



VIII. RESONANCES AT 0(e) 

The perturbation given in eq. (|22p shows that for each combination of n, to and k there result 8 harmonic vectors 
m. Again, due to the even character of the perturbation, we take just one representative resonant vector, m, whose 
coefficient a m also encompasses the contribution of its opposite vector — m; and keep only those harmonics such that 
0(a m ) < G(l/23 2 ). 

Since the harmonic vectors at 0(e, 1/23 2 ) can be splitted into two groups, we denote by y the subset of vectors 



10 

whose third component is zero and by Z that of vectors having their second component equal to zero. Therefore, the 
perturbation can be written in the fashion: 

V(I,0) = V 12 (I) J2 a m cos(m-6)+V 13 (I) ^ a m cos(m-6). (23) 

On applying the resonance condition, m ■ w = 0, with m £ Z 3 /{0}, to the unperturbed system, the following 
relation between the actions in each degree of freedom is obtained: 

mi ll /3 + m 2 ll /3 + m 3 ll /3 = 0, (24) 

so that each resonant vector could not have all its three components of the same sign. Furthermore, whenever a 
resonant vector has two of its components equal to zero, we get a null amplitude for the perturbation term. 
Thus, we obtain twelve different resonant vectors, grouped in the following set: 

V r (e,l/23 2 ) = {(2, -1,0), (2, -3,0), (2, -5,0), (4, -1,0), (4, -3,0), (6, -1,0), (6, -1,0), (2,0,-1), (2,0, -3), (2,0, -5), 
(4, 0,-1), (4, 0,-3), (6, 0,-1)} 

For those vectors m £ y, the resonance condition together with the energy conservation condition define a curve 
in action space (not just a point as in the 2D case) given by: 



[mi ] jr 

r„/ % ,13/4 (25) 

^ = {A-(i + ^) (/[)4 /3} , 



jr _ _ 



':! 



( 3/4 

+ (m 1 /m 2 )4 / 

The concomitant curve for m £ Z is given by 



where I[ £ [0,I^ iax ] with lf nax = U — — - 



v/A 



3/4 



*-{4-(l+3)«>-/ (26) 

rr _ _ f mi A jr 
J 3 - m 3 n' 



with I{ £ [0, /| m ], the upper bound being I^ lax = ( 1+(m i m u ) ■ In energy surface both kind of resonances define 



,3/4 

l+(mi/m 3 ) 4 / 

straight lines. On applying the pendulum approximation, we obtain a similar resonant Hamiltonian to that given by 
equation ((9]), namely 

H r (p 1 ,i> 1 ) = ^+eV m cosxl> 1 , (27) 

where V m = Vi2{I r )a m for m £y, and V m = Vi 3 ,(I r )a m for m £ Z. 

In energy space, the resonance widths in each degree of freedom, are adequately described by 

(Ah t ) r m = UitfXAIiYn = ^All /3 2(eMV m ) 1/2 mi = ^All /3 (eMV m )^ 2 m t . (28) 

Let us remark that, while in the 2D model the resonance and energy conservation conditions force the resonance 
width to depend on just one variable, either h, h\ or I[, in the 3D model the resonance width is a function of two 
independent variables, which we have chosen to be /{* and h. 

With the widths computed by means of (f2"8")l . we can trace the displacements of the resonant energies, h\ + (Aft,,)^ 
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and h\ — (A/ij)^, for values of I\ £ [0,I max \. Therefore, following [T(| we perform the global change of coordinates: 



ei = —=(h\- 2h 2 + h 3 ), 
V6 



(29) 



where ei € 



V 3"-' v^ 



e 2 € — -4=, -4= , e 3 = -4=, adopting the value h sa 0.485, to finally display in Fig. [6] 

the region of energy surface occupied by the structure of resonances at 0(e, 1/23 2 ) for two different values of the 
perturbative parameter. Let us remark that many of the resonances in V r (e, 1/23 2 ) are barely observable due to their 
thinncs and close proximity to a boundary. 



h=0.485, e=0.01 



h=0.485, e=0.05 




-0.4 -0.3 -0.2 -0.1 0.1 0.2 




FIG. 6: Resonances at 0(e, 1/23 2 ) for e = 0.01 (on the left) and e = 0.05 (on the right). 



IX. RESONANCES AT 0(e 2 ) 

As in the 2D case, we perform a canonical transformation in order to remove the perturbation at 0(e), the associated 
generatrix function being of the form: 



F(J,6) = J 6+ e$(J,0) 



where 



<S>(J,9) = J2 $™(J)sin(m-0)+ ^ $ m (J)sin(m • 0), 



(30) 



(31) 



with $ m (J) = -V m (J)/m ■ w(J). 

The new Hamiltonian is described by the same formal expression given in equation (| 1 3[) , where the first and second 
terms within braces adopt the values: 



l 6> 2 ff (J)d$(J,y)d$(J,y) 
2 dJjdJi dtpj dipi 



(32) 
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meym'ey 



I i ,2/3 ,2/3 . i ,4/3 x 

" ; T i/3 1 ; 7 i/3, i cos (( m + m ) ' y) + cos (( m -"»')•¥>)} + 



(miJ 1 + m 2 J 2 )(m' 1 Ji / ' 3 + m' 2 J 2 /0 ) 



E E **»■' 

E E 5m ' r 



/ , T 2/3 T 2/3 . , ,4/3x 

, 1/3 * '.iL ,,1/3' ,,1/3, l cos ((™ + ™ ) ' *>) + cos (( w - "» ) ' ¥>) } 
\m-\_Ji + mz J 3 ){m 1 J 1 +m 3 J 3 j 



1 7 2/3 7 l/3 7 l/3 

m 1 m 1 J 1 J 2 J 3 



mey m 1 eZ 



(mi J-/ + m 2 J 2 ' ^TOiJi + w 3 J 3 ) 



— jyj- {cos ((m + m/) • cp) + cos ((m - m/) • <p) } + 



__ / ,2/3 ,1/3 ,1/3 

+ EE <W' 7 — 7 i/ 3 mim 7 i/3w 2 , ,1/3 , — T7V3T i cos (( m + m ') ' ^) + cos (( m -"»')■ V) } . 

m££m'£)) (™lJl +^3^3 J^m 1 J 1 +m 2 J 2 ) 

where <5 m , m / = 2 3 3 2 / 3 /3 20/3 a rn ,ow and 
SV(J )V )9*(J,v) _ 



cu- 



efy?i 



(33) 



meym'ey 



25 r , 



E 

E E 2 * 

m-6.2 m' dZ 



(2m' 1 J 1 1/3 J 2 2/3 +m 2 J 1 4/3 J- 1/3 ) 



,1/3 



l/3x 



(2m' 1 J 1 1/3 J3 2/3 + m 3 J 1 4/3 J 3 - 1/3 ) 



,1/3 



(mur+™^ 3 /o ) 



l/3x 



{cos ((m + m/) • (p) + cos ((m — m') ■ ip) } 



{cos ((m + m') • ip) + cos ((m — m') ■ ipj } 



, ,1/3 ,1/3 ,1/3 
^ ^ 2 2 £ miTn / - -^3 1 2 ; 3 1/3 {cos ((m + m') • <p) + cos ((m - m') ■ v) } ' 



E E 2 ^ 



Therefore, the Hamiltonian may be recast as 



(miJi /3 + m 3 Jr) 
m > jV3 jl/3 jl/3 

■ l ' 2 - 3 1/3 {cos ((m + m') ■ y) + cos ((m - m') ■ tp) } . 



m ' m , , Tl/3 , , ,1/ox 

(m'lJi/ +m 2 J 2 ' ) 



H(J, ip) = H (J) + e 2 2_. \ E E [ a ( n ! m J m/ ! ^) + a(n,m, — m', J)] cos(n • <£>)+ 

neA I rreeCV rre'eCV 

+ 2_, /_, [^( n : m,m', J) + b(n,m, — m', J)] cos(n • ip)- 

me2. m' ez 

+ 2, /, [c(n,m,m' , J) + c(n,m, — m' , J)] cos(n • <p)- 



mey m' ez, 



\^ \^ {d(n,rn,rn',J)+d(n,rn,—rn',J)]cos(n-(p)>, 
mezm'ey ) 



(34) 



where A is the set of harmonic vectors arising through any combination of vectors from yuZ and whose first nonzero 
component is positive, and the coefficients a, b, c, and d are defined as follows: 



a(n,m, m', J) 



, I r 2/3 7 2/3 . / 7<l/3\ 

~l ,1/3 ,1/3 \ I / ,1/3 , ,1/3, 

(miJ 1 +m2J 2 )("li^i + rn 2"'2 ) 



-2 



I I ,1/3 , , ,1/3, 



if ± n = m + m/, 
if ± n 7^ m + m'. 



b(n,m, m', J) 



I 1 ,2/3 ,2/3 . / ,4/3n 

7 ,1/3 . ,1/3 w , ,1/3 . / ,l/3s 



(2rn[J^J^+rny^J£^) 

I I ,1/3 . / ,l/3x 

(m 1 J 1 > +m 3 J 3 ' ) 



if ± n = m + m', 

if ±n y^ m + m' . 
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c(n,m, m', J) = 



Jm.m' 



I ,2/3 ,1/3 ,1/3 

ynim 1 J 1 / J 2 ' J 3 

7 T l/3 ,1/3w , ,1/3 , 7~A 

(miJ 1 / +m 2 J 2 J^m^Jj +m 3 J 3 



, ,1/3 ,1/3 1/3 
I "'l^i '-*2 J 3 
7 / jl/3 i / ,1/3, 

(m 1 J 1 / +771373' ) 



if ± n = ra + m', 
if ± n 7^ m + m' . 



d(n, m,m',J) 



'm,m' 



1 ,2/3 ,1/3 ,1/3 

7ni777 1 J 1 / J 2 ' J 3 
(miJ 1 / +777.3J 3 ' )(777 / 1 J 1 / +777 2 J 2 ' 



-4 






7 / ,1/3 , ,1/3* 

KA + m 2 J 2 ) 



if ±n = m + m', 
if ± n ^ m + m'. 



Further, on introducing a coefficient anologous to the one in the perturbation at 0(e 2 ) for the 2D model, 
U)(n,J)=2, /_, [a{n,m,m', J) + a(n, m, — m' , J)] + y, /, [b(n,m, m' , J) + b(n,m, — m', J)] + 

+ 2. /, [c(n, rji,Tre', J) + c(n,m, —m', J)] + /, /, M( n i m,m', J) + d(n, m, m', J)] 

mey m'eZ m&Zm'ey 

the Hamiltonian allows for the expression: 



(35) 



H(J, tp) = H Q (J) + e 2 J2 P ( n ' J ) cos ( n ' f) 



neA 



(36) 



Close to any resonant action associated to a resonant vector at C(e 2 ), n, we could apply the pendulum approximation 
to get the same formal expression for the resonant Hamiltonian given by (|17j) . 
The variations of the energy components are similar to those of 0(e); indeed, 



Ah, 



Ujl /3 ^\Uo\) 1/2 ni 



For those n with either n^ — or ni = 0, the resonant curves are given by equations (j25|) and (|26|) respectively, 
with the pertinent substitution oinii and Ii by m and Ji. Meanwhile, whenever n\ = the resonant curves are given 
by 



J 2 (l+(n 2 /77 3 ) 4 ) 3/4 

J 3 = (l+("3/™2) 4 ) 3/4 •• ^ 



3/4 



{A_ (J[ )4/3 } 
{ ,_ (J[) 4/3 } 3/4 



(37) 



with jf e [o,(^) 3 / 4 ]. 

Let V r (e 2 ,l/23 2 ) be the set of resonant vectors which belong to A and that can be constructed by at least one 
pair (m,m') such that 0(a m a m ') < C(l/23 2 ). On computing the elements of V r (e 2 , 1/23 2 ) we learn that this set 
consists of fifteen vectors with one null component together with forty eight having all its three components different 
from zero. 

Fig. shows the area of the canonical energy surface (h — 0.485) occupied by those resonances at C(e 2 ,l/23 2 ) 
that have one null component together with the (2, — 1, —1) resonance, for the same two values of the perturbative 
parameter used for Fig. [6l 

Such a subset of V r (e 2 , 1/23 2 ) as well as the complete set V r (e, 1/23 2 ) have been considered in Fig. [U-left for 
e = 0.005. This picture should be compared with the contour-plot obtained by means of the MEGNO for the very 
same value of e, which is displayed in Fig. [8]-right (taken from [10j). This numerical exploration evinces that the 
resonances which strongly manifest are those with just one null component (i.e. the straight ones) and the (2, — 1, —1) 
resonance (the one showing a curved shape). 

A glance at Fig. [5]-left reveals that in some intersections between 0(e) and C(e 2 ) resonances, the widths of the 
latter tend asymptotically to infinity. This is due to the emergence of small denominators in the Fourier coefficients 
of the perturbation, a fact that reminds us that the canonical transformation perfomed in order to eliminate the 
perturbation terms proportional to e is no longer valid in the neighbourhood of any 0(e) resonance. 

A good example of this behaviour is the intersection between the (0,1,-1) and the (2,-1,0) resonances. The 
former is an 0(e 2 ) resonance that starts on the top right-hand corner of the energy surface and gets through the 
middle of it while the latter is of 0(e) and can be identified as the widest of the resonances departing from the bottom 
right-hand corner of the energy surface. 
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FIG. 7: Resonances at C(e 2 , 1/23 2 ) for e = 0.01 (on the left) and e = 0.05 (on the right). 
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FIG. 8: Resonances up to 0(e 2 , 1/23 2 ) for e = 0.005 (on the left) and the actual resonance structure obtained with the MEGNO 
(on the right) for the same energy normalized to h — 1. 

X. ANALYTICAL ESTIMATE OF THE CRITICAL VALUE OF THE PERTURBATION PARAMETER 

FOR THE 3D SYSTEM 



As it can be seen from Fig. [8l many resonances have a triangular shape. Such is the case of all those resonances 
associated to a vector having either its second or its third component equal to zero. 

It can be demonstrated instead that for a resonant vector with its first component n\ equal to zero, the coefficicntes 
a(n,m,±m',J) and b(n,m,±m',J) are null V (m,m'). Thus, the contribution of such a vector proceeds only 
through its concomitant coefficients c(n,m,±m',J) and d(n,m,±m',J). Further, from equation (|37|) it can be 
stated that the resonance width tends to zero when J\ approaches either or (h/A) 3 ^ 4 . Consequently, the regions 
encompassed by the separatrices of resonances for which m — do not have a triangular shape. 

On estimating e c for the 3D model, we are compelled to make somewhat strong simplifications: (i) we take as e c 
the value of e for which the total area covered by resonant regions (A r ) equals 90% of the whole area of the energy 
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surface (Ah); (ii) we approximate by triangles the resonant regions corresponding to resonant vectors with one null 
component; (hi) we approximate by two triangles the resonant region corresponding to the (2,-1,-1) resonance; 
(iv) we do not consider any further resonance; (v) we add up the area of each resonance disregarding the intersections 
due to crossings of resonances, so that those regions corresponding to two different resonances are considered twice. 
In Fig. [9]-left we have plotted the fraction A r /Ah for the perturbation parameter varying in the range e S 
[0.00001,0.2]. There it can be observed that A r (e) reaches 90% of A h for some e c between 0.03 and 0.04. This 
result is in quite good agreement with that arising from Fig. fright (taken from [12j) which displays the fraction of 
chaotic motion according to the MEGNO values, for the same range of e. 
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FIG. 9: 

scale 



Ar/Ah (on the left) and fraction of chaotic motion (on the right) both vs. the perturbation parameter, in logaritmic 



XI. DISCUSSION 



We have checked out the accuracy of the overlap criterion when applied to a simple near intcgrable Hamiltonian 
system in both its 2D and 3D version. To this end, we have computed the unperturbed resonances up to order 0(e 2 ) 
for both systems, and modelled each resonance by means of the pendulum approximation in order to estimate the 
theoretical critical value of the perturbation parameter for a global transition to chaos. 

By performing several surface of sections for the the 2D case we have derived an empirical value to be compared 
to our theoretical estimation, both being in good agreement. For the 3D case a theoretical estimate of the critical 
parameter has been attained, which is shown to match the one given in [12] . where such a value is achieved on 
computing the fraction of chaotic motion vs. e according to the MEGNO values. 

Let us remark that the conception of transition to global chaos assumed for the numerical estimate of e c in the 
2D case is of a different nature from the one adopted for the 3D system. Actually, the 2D system is considered to 
be globally chaotic if the chaotic component of phase space appears almost fully connected, while in the 3D case the 
system is regarded as globally chaotic when at most 10% of the energy surface corresponds to invariant tori. Notice 
that in the latter case, though it is very likely that the chaotic component be connected when resonances do overlap 
in a mostly chaotic phase space, nothing could be asserted about the existence of a fully connected region of unstable 
motion (see [12| , [3 and [l5| for a thorough discussion) . 

Therefore, from both theoretical and numerical results we may assert that a suitable estimate for the critical value 
of the perturbation parameter could be obtained by means of the overlap criterion when considering resonances up 
to 0(e 2 ). Indeed, regarding terms just up to 0(e) largely overestimates e c , as already shown by QJ for the Standard 
Map. 
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